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Abstract 

We investigate the totally asymmetric exclusion process on Z, with 
the jump rate at site i given by ri = 1 for i =/= 0, ro = r. It is easy 
to see that the maximal stationary current j(r) is nondecreasing in 
r and that j{r) = 1/4 for r > 1; it is a long outstanding problem to 
determine whether or not the critical value Vc of r such that j(r) = 1/4 
for r > Tc is strictly less than 1. Here we present a heuristic argument, 
based on the analysis of the first sixteen terms in a formal power series 
expansion of j{r) obtained from finite volume systems, that rc — 1 and 
that for r S 1, j{r) ~ 1/4 — 7exp[— a/(l — r)] with a sa 2. We also 
give some new exact results about this system; in particular we prove 
that j{r) — Jmax(?'), with Jmax('") the hydrodynamic maximal current 
defined by Seppalainen, and thus establish continuity of j(r). Finally 
we describe a related exactly solvable model, a semi-infinite system in 
which the site i = is always occupied. For that system, r^"' = 1/2 and 
the analogue j*'"' (r) oi j{r) satisfies j^"'(r) = r(l — r) for r < r^"'; j^"'(r) 
is the limit of finite volume currents inside the curve |r(l — r)\ — 1/4 
in the complex r plane and suggest that analogous behavior may hold 
for the original system. 



*Dedicated to Raghu Varadhan. This problem was part of a talk by JLL presented at 
the celebration of Raghu's seventieth birthday. We hope that this note will prompt Raghu 
to solve the problem. 
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1 Introduction 

Our understanding of nonequilibrium stationary states (NESS) of macro- 
scopic systems is very incomplete at present. While systems in equilibrium 
exhibit behavior which, away from phase transitions, is essentially local, the 
behavior in a NESS is global. That is, away from criticality, changes made in 
one region of an equilibrium system have almost no effect far away, but any 
local change in systems carrying fluxes of conserved quantities will almost 
certainly be felt at a distance. 

This difference in behavior is particularly striking in low dimensional 
systems, where local changes in conductivity can radically affect the total 
flux through the system. Some years ago one of us (JLL), together with 
S. Janowsky, studied a very simple system of this kind, a perturbation of 
the familiar one-dimensional totally asymmetric simple exclusion process 
(TASEP). Recall [HIS] that the latter is a model in which particles occupy 
some of the sites of a one dimensional lattice; we let r/j = 1 if a particle 
is present on site i, r]i = otherwise. The particles independently try to 
jump to the site on their right with rate one, succeeding if the target site is 
empty and otherwise staying at their original site. The lattice may be either 
(i) a ring, (ii) a finite interval with open boundaries which particles enter 
at rate a on the left and exit at rate /3 on the right, or (iii) all of Z. The 
invariant measures for these systems are known: (i) on a ring with L sites 
and particles all (^) configurations have equal probability; (ii) on the 
interval the system has a rich phase diagram in the a-f3 plane, obtained by 
the "matrix ansatz" [3] or otherwise [4]; (iii) on Z [T] there are translation 
invariant measures, the extremal ones of which are product measures Vp, 
< /O < 1, in which the probability of a site being occupied is p. In all 
cases there is a steady state current given by {r]i{l — ?/j+i)) (the value is 
independent of i), where for any function f{r]) of the configuration we write 
(/) for the expectation of / in the stationary state. The current in the state 

on Z is thus p{l — p), with maximum value 1/4 occurring for p = 1/2. 
(On Z there are also so-called "blocking" measures with zero current, which 
will play no role in our discussion.) 

Janowsky and Lebowitz [H [6] studied what happens in these systems 
under the introduction of a blockage: a change in the jump rate, from 1 
to some value r > 0, across one bond of the system. Consider first the 
system on Z, and suppose without loss of generality that the blockage is on 
the bond between site and site 1. We will be interested in the stationary 
states of the system and in particular in j{r), the maximum value, taken 
over all stationary states, that the current can assume. Clearly if r = then 
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j{r) = since all stationary states will have zero current: any particles in 
the system will just pile up behind and empty out in front of the blockage. 
For r > 1 the "blockage" in fact offers no impediment, in the sense that 
j{r) = 1/4 (see discussion below); more generally, one expects that in this 
case the blockage will have only a local effect and that, just as in the original 
(r = 1) system, there will be for any p £ [0, 1] a stationary state with current 
p{l — p) in which the measure far from the blockage approaches fp, and thus 
in particular that j{r) = 1/4. The question then is what is the nature of the 
stationary measures for r G (0, 1)? How will the system respond globally 
to this local change? It may be convenient to think in terms of traffic flow 
in which the width of one stretch of the road is changed. How well can the 
traffic adjust to a narrowing of the road? 

A simple mean field approximation gives a current r/{l + r)^ for r < 1, 
corresponding to Tc = 1 [6j. Beyond this, some precise information can 
be obtained by considering the system on Z as the infinite volume limit 
of the system either on the ring or on an interval. Consider in particular 
the open system on an interval of 2L sites labeled —L + 1, . . . , L, with the 
special bond between sites and 1 and with a,(3 > 1. It can be shown by a 
coupling argument [6] that the current j'^'^{r) in this system is decreasing 
in L for fixed r and, by a similar argument, that ji^'^{r) > j{r). A limiting 
NESS exists as L — )• oo (in theory one may need to take this limit along a 
subsequence, but that has no effect on the conclusions we draw below, and 
we ignore this possibility for notational convenience) so that 

j{r) = lim il'^ir), a,/3 > 1, (1.1) 

(no subsequence is needed here). Since clearly j'^'^{r) is nondecreasing in 
r at fixed L, j{r) must also be nondecreasing in r. Moreover, j(r) < 1/4 
for all r, as one may see by an argument based on [2l p. 273]: consider 
an unbounded sequence of translates of the stationary state for which a 
limiting state exists; this will be a steady state for the unperturbed system 
on Z with all rates equal to 1 and thus have current at most 1/4. Since 
= 1/4, there exists a critical value Tc S [0, 1] with j(r) < 1/4 for r < rc 
and j(r) = 1/4 for r > rc- The key question that we would like to answer is 
whether or not Tc = 1. 

Exact formulas for ji'"{r) for L = 1, 2, 3 and for j^'^{r) were obtained in 
[6]. These show that j4'^(0.51) < 1/4, which implies that rc > 0.51: this is 
unfortunately the best rigorous lower bound that we know, although numer- 
ical simulations strongly support the conclusion that Vc > 0.75. Simulations 
on large open systems in [6] gave a smooth concave curve for ji^{r) ~ j{r), 
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r G [0,1] which appeared to approach the value 1/4 very flatly. Janowsky 
and Lebowitz also proposed an analytic form which fits this data very well: 

i(r)« 1/4 -exp[-G(r)/(l-r)/7(r)], rE [0,1], (1.2) 

with G and H polynomials determined (as Pade approximants) from the six 
terms in the power series J(r) given in (j3.ip below. 

There has been further rigorous work on this and related problems. 
Covert and Rezakhanlou [7] study the hydrodynamic limit of a system in 
which the site-dependent jump rate is a discretization of a uniformly positive, 
bounded, and continuously differentiable function r(x), x G M, and prove 
the existence of a hydrodynamic scaling limit in which the evolution of the 
density p(x, t) is governed by the Burgers' equation dtp + dx[rp{l — p)] = 0. 
Seppalainen [8] proves the existence of the hydrodynamic limit of the system 
with a single slow bond, as considered here. He shows also that the maxi- 
mum current Jmax('^) in the model, as he defines it, is a continuous function 
of r, but leaves open the question of whether or not his and our defini- 
tions of the maximal current agree. We sketch in Appendix |A] an argument 
establishing that Jmax(?') = jif), so that in particular j(r) is continuous. 

As to the value of Tc, it is argued in [11], from a finite size scaling 
analysis of simulation data, that rc ~ 0.80. On the other hand, our process 
is directly related to a certain last passage percolation model [Hl[9l[T0], and 
this in turn is equivalent to a directed polymer model on the lattice, at 
zero temperature, with a line of defects, that is, with a columnar "pinning" 
potential; see |iOj for a review. In this language Tc = 1 corresponds to 
the pinning of the polymer for arbitrarily small values of the potential. 
The current interest in such systems, which have been studied in both the 
physics and mathematics literature, is one motivation for the study of the 
blockage model. In the physics literature, renormalization group arguments 
[T2j indicate that a directed polymer at low temperature in the continuum is 
indeed localized for all values of the potential, giving some support to Tc = 1. 
Mathematically, Beffara et al. jTSj study a model of random polynuclear 
growth with a columnar defect and find a phase transition at a nonzero 
value of the pinning potential, although one of the authors believes that 
Tc = 1 for the blockage problem |14j . 

The present paper is primarily concerned with the value of rc and the 
nature of the limit j{r) — )• 1/4 as r — t- Tc. We give arguments that Tc = 1 
and that for r ^ 1, j(r) ~ 1/4 — aexp[—A/{l — r)] with j4 ~ 2, a form 
consistent with (jl.2p . These arguments are far from a proof; they are based 
on the analysis of sixteen terms of a formal power series whose convergence 
is not known but whose derivation suggests that it may represent j(r). 
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Remark: The system on a ring. For computations it is frequently 
convenient to consider a finite system on a ring rather than on an interval; 
we always consider a ring with 2L sites, labeled as for the interval above 
and again with the special bond joining and 1, and with L particles. We 
believe on the basis of the behavior when r = 1, the exact solution of small 
systems, and the simulation of large systems, that the current j"°^(r) in this 
system is decreasing in L and satisfies limi_!.oo = but we do not 

have a proof. A coupling argument leads to the inequality j^'^ {r) > j""^(r) 
for a, /3 > 1. 

The rest of the paper is organized as follows. The derivation of the power 
series referred to above is presented in Section [2] and its analysis is given in 
Section [3l In Section [4] we discuss a simpler model, a semi-infinite TASEP 
on sites 1,2,... in which particles enter site 1 at rate r; here the relation 
between the current and its power series can be completely determined, and 
we argue that the results shed some light on the properties of the original 
system. Some technical points are relegated to appendices. 



2 Power series 



The steady state probabilities in any of the finite systems considered above — 
on the ring or the interval — are rational functions of r, since they may be 
obtained as ratios of minors of the generator M(r), which is affine in r (see 
()2.6p below). The current is thus also rational in r; it is nonsingular (and in 
fact vanishing) at r = 0, and so has a power series in r convergent in some 
neighborhood of the origin: 
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3l'\r) = Y: cli r^ (r) = ^ vK (2.1) 

fc=0 fc=0 

It is shown in |6j for the system on the interval, and may be verified similarly 
for the half-filled ring, that for k < L the coefficients in ()2.ip agree and are 
independent of L, a, and /3 as long as a, /3 > 1; we denote the common value 
as Ck- 

Ck = cti = clJ, k<L, a,/3>l. (2.2) 
Thus for jl = jl'^ or H = j™^, 

L 

jl{r) = Y,cur^ + 0{r^+^), (2.3) 

k=Q 
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where the remainder term, but not the c^, may depend on L and on which 
current is under consideration. 

More generahy, one can by similar methods obtain a corresponding result 
for the entire measure. The proof of the next result, and of the extension 
mentioned in Remark 12.2( a) , is implicit in the proof of Theorem 1 in Ap- 
pendix B of [6] . 

Theorem 2.1. Suppose that L > j > 1, and let f{r]) be a function of the 
configuration which depends only on the 2j variables ?7_j+i, . . .r]j. Then in 
each of the finite systems we are considering (/) has a power series in r for 
which the coefficients up to order L — j are independent of L. Moreover, on 
the interval these coefficients are independent of a and f3, for a,l3 > 1, and 
are the same as on the ring. 

Remark 2.2. (a) In the special case of an open system with a = (3 = 1 one 
in fact obtains one more order, that is, the coefficients of the power series 
of (/) up to order L — j + 1 are independent of L. 

(b) (|2.2p is a corollary of the j = 1 case of Theorem 12. H since the current 
in any of these systems has value r (r/o (1 — r/i)). 

We may compute the power series (j2.ip . or indeed the series for any 
(/) on the ring or interval, by a direct expansion, without obtaining the 
exact rational form (which is usually a much more computationally in- 
tensive problem). For if the steady state measure assigns probability 
P{v) = Yl'k=oPk{'n)'''^ to the configuration r] then P is the null vector of 
the generator M{r) of the process: 




(2.4) 



M(r) has the form M(r) = M^^^ +rM^^^ and po{rj) is known, so that pi{r]), 
P2{'n)-, etc. may be found recursively by solving 



M^p;- = -M(Sfc_i. (2.5) 



For example, for the ring po(^) = ^nr)(°) with r/^ = 1 if i < and = 
ii i> 1, and M is a i^j^^ x {^j^^ matrix with 



'77,(1 -7?,+i), if C = 7?^'^+\ i ^ 0, 
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Here rj^^ denotes the configuration obtained from r] by exchanging the values 
of r]i and r/j. M^^^ is not invertible but ()2.6p implies that it has a triangular 
structure so that (|2.5p may be solved by back substitution. We have used 
this procedure to determine the coefficients of (j2.2p to order 16 and of 
the densities (r/j), i = 1,...,5, to order 16 — i; the results are given in 
Appendix |BJ The densities behind the blockage may be obtained from this 
table and the symmetry {r]i) = 1 — (r/i_j). Although the Ck given in the 
table alternate in sign we believe, on the basis of an argument given below, 
that this alternation will break down at A; = 17. 

In the next section we make, for the purposes of discussion. 

Assumption A: The series YlT=o '^^f^ ^ nonzero radius of convergence 
and represents the current j{r) for those r > for which it converges. 

We ask what this assumption, together with the knowledge of the first six- 
teen Cfc, suggests about the behavior of j{r). Such a procedure is clearly 
fraught with uncertainty, and any conclusions drawn must be approached 
with skepticism, both because the series might diverge or, if convergent, 
might represent j{r) only in a limited region or not at all, and because 
knowledge of any finite number of coefficients of a power series can tell us 
nothing with surety of the behavior of the sum. Nevertheless we feel that 
there is justification for our approach: the analysis of the series seems fairly 
compelling, and the conclusions we reach are consistent with the findings of 

3 Analysis of the series for the current 

Let us define 

oo n 

J{T) = Y,Ckr\ Jn{r) = Y.Ckr\ (3.1) 

fc=0 fc=0 

where the coefficients Ck are defined by (12. 2p and cq , . . . , cie are given in 
Table [D of Appendix [HI as indicated above, we assume for the purpose of 
this discussion that the series for J(r) has a positive radius of convergence. 
We let J(r) also denote the analytic continuation of the sum of this series 
to other values of r. The key tool in our analysis will be the fact that for 
any function f[z) analytic at z = the Taylor series of F{r) := f{J{r)) in 
powers of r can be calculated up to order 16 from the known coefficients c^; 
we will write -Fi6('^) for the resulting polynomial. 

We first note that for k = 1, . . . , 16 the coefficients alternate in sign 
and that for values of k in the range 3 < A; < 13 the ratio Ck/ck-i varies in 
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magnitude in a reasonably narrow range around 2/3. We accept the natural 
explanation of this behavior: that J(r) has a simple pole with relatively large 
residue at a point r = rg ~ — 1.5, affecting the early coefficients of the series, 
but that this effect is masked in later coefficients by weaker singularities 
closer to the origin. Further analysis of the series will be aided by a more 
precise location of tq. There are various methods to achieve this; we discuss 
two of them in AppendixO The conclusion is that we may take ro = —1.5437 
and expect this value to be accurate within an error of ±0.0005. In the 
remainder of this section we will often study the power series for (r — ro) J(r) 
(or some similarly constructed function — see, e.g., (j3.3p ) rather than that 
for J{r), in order to eliminate the strong effect of this pole and thus make 
other singularities of J(r) more apparent in the behavior of the power series. 

We now consider the question of whether or not rc = 1 . Under Assump- 
tion A, the radius of convergence of J(r) can be at most rc, so that the 
Taylor coefficients of J{r), or of (r — tq) J(r), must grow (for large k) at 
least as fast as r~'^. Figure [1] shows the latter coefficients for = 6, . . . , 16; 
there is no sign of exponential growth, suggesting that rc = 1. 



0.02- 



0.01- 



1 1 1 — 

10 12 



14 16 



-0.01- 



Figure 1: The Taylor coefficients of (r — ro) J(r). 

Accepting provisionally that rc = 1 we now ask about the form of the 
singularity in J(r) there. The heuristic formula (|1.2p suggests that 

J(r) ~ 1/4-7 e~"/(^-'') for r < 1, 7, a > 0, (3.2) 

so that (1/4 — J{r))^^ ~ 7~^exp(a/(l — r)). To verify the consistency 
of this ansatz we observe (see Appendix |D|) that the Taylor coefficients of 
exp(o/(l — r)) grow as const • fc-3/4g2v^(^]^ _j_ o(l)). With Uk the Taylor 
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coefficients of (1/4 — J{r))^^ we use the function Statistics[Fit] in Maple to 
obtain a fit (also shown in Figure [2]) over the range 7 < /c < 16 of the form 
Uk ~ AiVk + Bi log A; + Ci with Ai = 2.82, 5i = -0.495, and Ci = 0.116; 
the residual standard deviation in the fit, that is, the square root of the 
ratio of the residual sum of squares to the number of degrees of freedom 
(here 7), is 0.0002, quite small. The value of Ai is consistent with the value 
o ~ 2 which we discuss below; Bl, CI (and higher terms if sought) are 
more sensitive to corrections to the exact type of the singularity, to other 
singularities, etc. 

100- 
90- 
80- 
70- 
60- 
50- 
40- 
30- 
20- 
10- 

2 4 6 8 10 12 14 16 

Figure 2: A plot of log^ Uk together with a straight line fit, where the Uk are 
the Taylor coefficients of (1/4 — J(r))~^. 

To explore (13. 2p further we consider 

X(r) =log[(r-ro)(l/4- J(r))]. (3.3) 

If ()3.2p holds then X will have a simple pole at r = 1, which by itself would 
lead to a Taylor series for X with all coefficients equal to —a; due to the 
presence of other singularities we would of course expect such constancy at 
best asymptotically. A singularity of X for |r| < 1 would lead to exponen- 
tially growing coefficients, while other singularities on the unit circle would 
modify the asymptotic form. 

The Taylor coefficients Xk of X are plotted in Figured There is no sign 
of exponential growth, again suggesting that Tc = 1 (and more generally 
that J(r) 7^ 1/4 for |r| < 1). At a first approximation one has Xk ~ —2, 
which is consistent with (|3.2p with a ~ 2. Overall, we believe that this 
analysis gives good support for the conclusion that Tc = 1 and that j(r) has 
a singularity of the general form (13. 2p at r = Tc = 1. 
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Somewhat more speculatively, one may in fact obtain a more detailed 
hypothesis for the singularity structure of J(r). Figure [3] suggests an os- 
cillatory behavior rather than constancy for the x^- Using the function 
Statistics[Fit] in Maple we obtain a very good fit (also shown in Figure [3|) 
to a function of the form x{k) = A + Bcos{C{k — D)), with 

A = -2.00922, B = 0.193059, C = 0.260931, D = 0.919233, (3.4) 

over the range 7 < A; < 16. The residual standard deviation is 0.00025, 
about one hundredth of one percent of the typical value being fitted. 

k 




Figure 3: The Taylor coefficients of log[l/4 — (r — ro)J(r)], together with 
the curve A + B cos{C{k — D)) with parameters given by (j3.4p . 



The function X{r) with Taylor series given exactly by X]fcLo'^(^)'' 

A B[cosCD-rcos{C{D + m 
^ ^"^^ " 1^ + l-2rcosC + r2 ' ^"^'^^ 

with simple poles at r = 1 and r = e^**^. This suggests that one may well- 
approximate X{r) by the function Y(r) whose Taylor series agrees asymp- 
totically with that of X{r) but has the same first 16 terms as does that of 
X{r), 

16 

X{r) « Y{r) := X{r) + ^(xfc - x(A:))r^ (3.6) 

k=0 

and correspondingly approximate 
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By construction the first 16 coefficients in the power series for K{r) agree 
with those of J{r); the 17^^ coefficient is negative, suggesting that the al- 
ternation of sign noted in Table [1] breaks down at this point. 

One must again treat this analysis with caution; certainly for example 
the exact values of A, . . . , D depend on the range of indices k over which 
the fit is performed. 

Remark 3.1. One can attempt a similar analysis of the power series for 
the densities Pi{r) for small values of i, whose coefficients are also given 
in Appendix B. It is natural to guess that, near r = 1, pi would have a 
decomposition similar to (|3.2j) : pi{r) = pl^^{r) + /9^™^(r), with p^'^^ regular 
at r = 1 and pf^^ having an exponential singularity e~^^^^~^\ (Since the 
current is given by j{r) = r(r/o(l — = (??i(l — ^ / 0, where the 

expectation is taken in any NESS with maximal current, occurrence of the 
exponential singularity of (j3.2p in the current implies that it must occur 
in densities and/or in the two-body correlation functions.) The key to the 
analysis of the current, however, was that we knew (or could hypothesize) 
that the regular part of the current at r = 1 was precisely 1/4, and could 
therefore isolate the singular part for further study through the definition 
p.3p of X{r). On the other hand, while we know that Pi{l) = p^^^{l) = 0.5, 
it appears that the derivatives of pi (or equivalently of pj'^^) at r = 1 are 
nonzero, and we have not been able either to determine reliably any terms 
beyond the constant one in the expansion of pY'^ around r = 1 or to isolate 
pf^^ so that its behavior could be studied. 

4 A solvable semi- infinite model 

We would like to understand the convergence of the finite volume current 
j™^(r) (or j'l^'^ {r)) as L — t- oo, with particular attention to the behavior in 
the complex r plane, and for this purpose consider a semi-infinite system 
in which this behavior can be determined completely, using the work of 
Derrida, Evans, Hakim, and Pasquier [3j. To obtain this model we begin 
with the open system as described in Section [U with (3 = 1, and then replace 
the portion of this system lying to the left of the distinguished bond by a 
reservoir of particles, so that the site immediately to the left of this bond is 
always full. Now we study the portion of the system to the right of the bond; 
particles enter this system from the left at rate r and (in the finite system, 
now of L sites) exit on the right at rate 1. The current in this system is an 
upper bound for that in the original model [HI [2]. 
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The system may alternatively be viewed as a symmetrized TASEP on Z 
|10j : the configurations are assumed to be symmetric under the joint action 
of the reflection i — >■ 1 — i exchange of particles and holes, and the dynamics 
are such that a transitions occur simultaneously on symmetrically located 
bonds, at rate r on the bond from to 1, which is invariant under the 
symmetry, and at rate 1 on all others. Note that if the symmetry condition 
is satisfied by the initial condition then it is preserved by the dynamics. In 
this language the model is again related to a last passage percolation model, 
a discrete-time version of which was studied in |15j . 

The current j^^ in this system is computed in [3]: 

jVir) = ^1^, (4.1) 
where the polynomial Q^' has order L and is given explicitly by 

and implicitly through the recursion 

QlUr) = {1 - r)Qr{r) + j^(^^l^y^+\ L > 1, (4.3) 

with Qo'{r) = 1. It follows from (gSD that 

j|-i(r) =r(l-r) + 0(r-^+2). (4.4) 
It is further shown in [3j that 

/»=lim,r(r)=^(^-^)' if0<.<l/2, 

^ ^ L^oo-^^ ^ ^ |l/4, if r > 1/2. ^ ^ 

Thus the critical value of r in this model, r^"' = 1/2, is very different from 
Vc, and moreover the behavior of j^~^{r) near r^"' is very different from that 
conjectured in (13. 2p for j{r) near Vc- 

The analysis which leads to ()4.5p can be carried out also in the complex 
r plane. The simple closed curve T defined by |r(l — r)\ = 1/4, Re r < 1/2, 
divides the plane into an interior region ili and an exterior region Q2] see 
Figure m Then for complex r, 

•si/\ 1- -si/x 1^(1-''); lire 0.1, ,^ 
f\r) = hm jl\r) = { " ^ ' (4.6) 

L^oo 1/4, if r G 0,2- 
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We see also that is continuous on T. 

The mechanism for the convergence ()4.6p can be explored since the zeros 
of the polynomial Q^£^ can be obtained and plotted numerically for reason- 
ably large systems; this "Lee- Yang" approach to non-equilibrium systems 
was introduced in [16] and has been applied to the open TASEP in fL7\ . 
The zeros appear to fall on a smooth curve which approaches F as L in- 
creases; the limiting curve T and the density of the zeros were computed in 
|17j . The distance of the zeros from T decreases approximately as 1/L, and 
the rightmost (complex conjugate) zeros lie at a distance of order approxi- 
mately from the point r = 1/2. See Figured! which plots the curve F 
and the zeros of Q^' for L = 5, 10, 20, 40, and 80. These zeros are of course 
poles of j'l^"'. Note from (j4.4p that there is no trace of these singularities in 
the first L + 1 terms of the power series for j^^; rather they are hidden in 
the 0(r^"^^) remainder. 

Comparison of the models: We now want to describe a certain paral- 
lelism between the blockage and semi- infinite models, and to suggest that 
this may extend to some properties which can be established for the lat- 
ter but not for the former. To make specific comparisons, we will consider 
the finite-volume version of the blockage model to be that on the ring, but 
similar conclusions hold for the interval. 

Before discussing this parallelism, however, we note a key difference be- 
tween the models: in the representation ()4.ip of the numerator and 
denominator are polynomials in r of degree L, while if for the finite- volume 
blockage model we write 

ring -Pl^^C^) 

Ql i^) 

with P^™^(r) and polynomials, then the degree of P^™^ and Q^"^^ 

grows exponentially with L. This of course makes it difficult to carry out 
computations in the complex r plane for the blockage model similar to those 
described above for the semi-infinite model. To illustrate the difficulty we 
have calculated Q^l^^{r) for L < 5; the degrees of these polynomials are 1, 
2, 5, 14, and 42 for L = 1, 2, 3, 4, and 5, respectively. Their zeros are shown 
in Figure El and one may compare the complexity of this figure with the 
relative simplicity of Figure [H 

We next note two particular similarities between the models. 

1. In each model the coefficient of order k of the Taylor series for the 
finite volume current, fl^^{r) or jf^^r), is independent of L for L > k; 
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Figure 4: Zeros of the denominator Q^£^{r): The symbols are: L = 5, red 
cross, L = 10, magenta circle; L = 20; blue diamond; L = 40, green square; 
L = 80, black star. The solid line is the curve F = {r | |r(l — r)| = 1/4}. 



compare (j2.3p and ()4.4p . The parallelism is not complete: if we pass to 
the power series obtained from the L-independent terms we obtain in 
the semi-infinite model just a polynomial of degree two, r(l — r), and 
in the original model the rather complicated function J(r) studied in 
Section [2l but for our present purposes we will ignore this distinction. 

2. Although we have emphasized above that the singularity structure of 
j™^(r) is much more complicated than that of there neverthe- 

less are similarities. To see this, consider the two plots of Figure El 
Here we show: 



• for the semi-infinite model (Figure (D^a)) the zeros of Ql' for L = 
1, 2, 3, 4, 5, as well as the limiting curve F; 
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Figure 5: Zeros of the denominator (5^°^(r): L = 1, magenta circle; L = 2; 
blue diamond; L = 3, green square; L = A, black star; L = 5, red cross. An 
expanded view of the dashed square is given in Figure [6jb) . 

• for the blockage model (Figure EJb)), only the zeros of Q^l^^, 
L = 1, 2, 3, 4, 5, which lie near the origin; for the window chosen, 
exactly L zeros appear for (as for the semi-infinite model) 

except for L = 5, for which two extra zeros, farthest from the ori- 
gin, also appear. Also shown is a possible limiting curve discussed 
below. 



The similarity of these sets of zeros is clear. 

This parallelism suggests the following possible behavior for the blockage 
model: 

• Some subset of the zeros of Q^"^^ converge as L cxd to a curve T 
somewhat similar to T. The simplest hypothesis is that F is determined 
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Figure [6]^ a) 




Figure Mjo) 



Figure 6: Zeros of the denominator of the current for the semi- infinite model 
(a) and the blockage system on the ring (b), for = 1, magenta circle; 
N = 2] blue diamond; = 3, green square; = 4, black star; N = 5, red 
cross. Also shown are the limiting curve F for the semi-infinite model and 
an approximation to the conjectured limiting curve F of 
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by the relation 

\J{r)\ = 1 (4.8) 

in parallel with the definition of T. The curve shown in Figure [6|^b) 
was obtained by numerical solution of K(r) = 1/4, where K(r) is the 
approximation to J(r) defined in (|3.7p . 

• In the interior of the curve F, J(r) is analytic and the finite volume 
current j^l^^{r) converges to J(r) as L oo. In particular, Assump- 
tion A of Section [2] is valid. 

• The behavior of J(r) outside F has no direct connection with the 
behavior of the finite volume current there. 

This suggested behavior is not sufficiently well founded to be called a con- 
jecture, but we hope that it may be a guide to further investigation. 
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A Agreement of Jmax(^) and j(r) 

In this appendix we sketch a proof that the current j(r) defined in Section [1] 
is identical to the current Jmnxi'r) introduced in |8| and shown there to be 
continuous, showing that j{r) is continuous. 

Consider three versions of the 1-d TASEP process: r]t is a process on 
Z with initial configuration r/o(i) = 1 if i < 0, 7?o(^) = if i > 0; is 
likewise a process on Z but with arbitrary initial configuration (q, and is 
a process on the interval {—L + 1, . . . ,L} with entry and exit rates 1 and 
with initial configuration the restriction of r]Q to this interval. Transitions in 
these processes occur at rate 1 on all bonds except {0, 1}, where the rate is r; 
moreover, the dynamics are coupled so that the same Poisson clock on each 
bond controls transitions in each process, and the clocks on {—L, —L + 1} 
and {L,L + 1} control also the entry and exit, respectively, of particles in 
the T process. 
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Suppose now that /i is a measure on {0, 1}^ which is invariant for the 
TASEP dynamics with the rates given above and has (maximal) current 



j{r) = r(r7o(l — ^i))^- (^^o^ the sake of precision we define also j^^ir) = 
r(?7o(l — with fi^£^ the invariant measure for the open system.) Let 

Nj^{t) be the number of particles that cross the special bond in the rj process 
between times and t, with -/V^(t) and Nr{t) defined similarly. From the 
coupling one sees easily that Nrit) > N^{t) > N,^{t) for all t and all Coi so 
that for t > 0, 



E 





> E 




> E 




t 




t 




t 



(A.l) 



where E[-] denotes the average over the dynamics and in (-/V^ (t))/^ the initial 
configuration (^o is averaged over fi. But also 



•1.1^ \ 



E 



lim E 

t—>-oo 



Nrit) 



t 



t > 0, (A.2) 
(A.3) 

,(r); with (fTT]) and ([A?T]) - 



We show below that limt_^oo lE[A^,,(t)/i] 
()A.3p this completes the proof. 

For n > 1 let be the time of the n^^ crossing of the {0, 1} bond by a 
particle in the r]t process. It is shown in [8] that K{r) = lim„_j.oo Tn/n exists 
almost surely and is constant. The current Jmax is defined to be 1/K(r), so 
that Jraax{r) = linin^oo n/Tn a.s. It is easy to see that Nrj{t) — )• oo and so 



'Nr,{t) 



OO a.s. as t — >■ oo; from this and the inequalities 

Nr,it) + 1 



(A.4) 



we see that lim(_j.oo ^ri{i)/t = Jma.x a-S- We complete the proof by verifying 
that the family of random variables {Ny^{t)/t | t > 1} is uniformly integrable, 
that is, that 

■N,{t) 



lim max E 

C/-00 t>l 



t 



(A.5) 



where xs denotes the characteristic function of the set S. 

Note first that T„ > Sn, where Sn is the sum of n independent exponen- 
tial random variables of mean 1 (these can be taken, for example, to be the 

of [8j). Then for c > 2, and using ?7v^(t) < t, we 



variables Ya, 1 < i < n 
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have 



E 



Nr,{t) 



X{N^it)>ct} 



oo 

EM 



n 



n=2 

oo 

< 

n=2 



X{Nr,{t)=n}X{n>ct} 



n 



'X{n>cSn} 



(A.6) 



Now n/^n is gamma-distributed with density n^x" '^^/(n — 1)!, so that, 
using n! > n^e"", we have 



E 



n 



"X{n>cS'n} 



< 



n 



n /"l/c 



(n-1)! 70 



x"- dx < 



(A.7) 



and (|A?5]) follows from ([XG]) and (A7i}i . 



B Tables of coefficients 

We give in Table [1] below the coefficients in the power series (j3.ip for the 
current J{r), to order 16, and in Table [2]the coefficients in the series 

oo 

Pi{r) = ^dikr'' (B.l) 

k=0 

for Pi, up to order 16 — i, for i = 1, . . . ,5. In Table [J we report the values 
to 20 decimal places. We believe that these digits are reliable because in 
our computation of cq, . . . C15 the first 24 digits obtained were stable under 
increase in the system size (this test is not available for ciq). For cq, . . . C13 a 
check for stability under change in the precision used in the program found 
the same reliability. 

C The value of rg 

As observed in Section [3l the power series for the current suggests the ex- 
istence of a pole at some value tq ~ —1.5. In this appendix we attempt to 
determine an accurate value for tq by two independent methods. 

Method 1. We let W{r) = (r — tq) J(r) and compute, for various values 
of ro, the Taylor polynomial WiQ^r) = X^^Lo '"^fc''^- '^^^ oscillation in the 
coefficients Wk should be minimized for the correct value of ro; to emphasize 
the oscillation we compute the second difference Wk = Wk+2 — 21^^+1 + as 
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k 


Cfc 




k 


Cfc 





0.00000000000000000000 




9 


0.09696668201649961665 


1 


1 .00000000000000000000 




10 


-0.05607405058503243547 


2 


-1.50000000000000000000 




11 


0.04033294103506195933 


3 


1.18750000000000000000 




12 


-0.02482314806701423240 


4 


-0.77889901620370370370 




13 


0.01477809455732788252 


5 


0.52961027553406380931 




14 


-0.01328263357536883488 


6 


-0.32787247554422253211 




15 


0.00277198065020739725 


7 


0.22700745336484005616 




16 


-0.00936905520626202337 


8 


-0.13514784111152134747 









Table 1: The coefficients of the power series J{r). 



k 


dik 


d2k 


dsk 


dik 


d5k 





0.00000000 


0.00000000 


0.00000000 


0.00000000 





00000000 


1 


1.00000000 


1.00000000 


1.00000000 


1.00000000 


1 


00000000 


2 


-0.75000000 


-0.68750000 


-0.65625000 


-0.63671875 


-0 


62304688 


3 


0.45312500 


0.40075231 


0.37939743 


0.36763019 





36002612 


4 


-0.32345016 


-0.31775049 


-0.32035731 


-0.32441081 


-0 


32884964 


5 


0.19113754 


0.16723143 


0.15934419 


0.15807536 





16042006 


6 


-0.13508488 


-0.13527025 


-0.14002291 


-0.14544675 


-0 


15090756 


7 


0.08218276 


0.07251374 


0.06758191 


0.06440565 





06221118 


8 


-0.05412219 


-0.05222456 


-0.05418729 


-0.05745286 


-0 


06108373 


9 


0.03743811 


0.03658774 


0.03602396 


0.03505521 





03382992 


10 


-0.01961547 


-0.01492080 


-0.01314674 


-0.01296838 


-0 


01376695 


11 


0.01873543 


0.02252206 


0.02556734 


0.02764994 





02887908 


12 


-0.00544623 


0.00051738 


0.00480082 


0.00795647 






13 


0.01025427 


0.01536259 


0.02003591 








14 


-0.00040069 


0.00489313 










15 


0.00558944 













Table 2: The coefficients dik of the power series pi = J2k dik^^ ■ 



well as w{k), the least-squares linear fit to Wk over the range A; = 8 . . . 14, and 
plot in Figure[7]the difference Wk — w{k) for ro = —1.5435 and —1.5440. The 
oscillations clearly change sign over this range of rg and will be minimized 
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at some intermediate value. 



0.00005 



-0.00005 



ft 



Figure 7: The difference Wk — w{k) for tq = 3.5435 (red squares) and rg = 
3.5440 (blue diamonds). 

Method 2. Consider the function V{r) = 1/(1/4 — J(r)), which has a 
zero at r = ro and for which we expect a singularity at r = 1, due both 
to the vanishing of the denominator and to the (hypothetical) exponential 
singularity of J(r) suggested in <\1.2\ . The fractional linear transformation 

u = u{r)=[- —]-^, ^ r = r{u) = -, (C.l) 

\ ri / 1 — r 1 + ri[l — u) 

with ri = —1.5, leaves the origin unchanged, carries r = ri to n = 1, and 
moves the singularity at r = 1 as far away as possible, to u = oo. We 
write V{u) = V{r{u)) and find that Viq{u) has one zero uq very near to 
u = 1, corresponding to ro = r{uo) = 1.5437, and that all other zeros 
are a distance at least 1.39 from n = 0. The value of ro found by this 
procedure is independent, to very high precision, of the choice of ri in (jC.ip 
for —3 < ri < — 1. 

Thus the two methods above give good agreement and suggest that ro G 
[—1.540, —1.535]. We take ro = —1.5437 and expect this value to be accurate 
within an error of ±0.0005. 
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D Asymptotics of Taylor coefficients 

In this section we find the large k asymptotic behavior of the coefficients 
in the Maclaurin series 

oo 

e»/M=^6fer^ |r|<l, (D.l) 

k=0 

where a is a positive number. We have 

oa/(l-s) 



h 



L / ^^^ds= I e-'^/*-('=+^)'°(i+*)di, (D.2) 



with C\ a circle of radius R <1 centered at 0, traversed anticlockwise, and 
C2 = -1 + Ci. The saddle points of -a/t - {k + \) ln(l + t) are: 



a + v^a(a + 4fc + 4) a- ^ a{a + Ak + ^) 
t+ = — T, ^ and t- = — . (D-3) 

The steepest variation path homotopic to C2 and passing through t- is an 
analytic curve orthogonal at t_ to the real line. The saddle point method 
applied at f = f_ gives 



6fc = const • fe-3/4e2^/«'=(l + o(l)) as A; ^ 00. (D.4) 
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